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FOREWORD 


The  experimental  data  originating  in  this  report  were  obtained 
by  J.  H.  Panesci  and  R.  C.  German.  Theoretical  results  were 
obtained  from  an  IBM  7070  computer  programmed  by  D.  H.  Hoyle. 


ABSTRACT 


The  two-dimensional  supersonic  base  pressure  theory  developed 
by  Dr.  H.  H.  Korst,  which  applies  to  two-dimensional  systems  with 
straight  jet  boundaries,  has  been  modified  to  be  more  applicable  to 
axisymmetric  ejector  systems.  The  modification  consists  of  a  new 
theory  for  estimating  the  peak  recompression  static  pressure  which 
is  applicable  to  either  axisymmetric  or  two-dimensional  jets  having 
straight  or  curved  boundaries.  Significantly,  the  recompression 
mechanism  is  not  independent  of  viscous  effects  for  systems  which  pro 
duce  a  non-uniform  inviscid  flow  field. 

Both  the  recompression  theory  and  the  modified  base  pressure 
theory  are  experimentally  verified  for  the  case  of  isoenergetic  mixing 
and  negligible  initial  boundary  layer.  The  recompression  theory  is 
shown  to  agree  with  experimental  results  to  within  ±10  percent,  where¬ 
as  the  experimental  base  pressure  results  have  a  standard  deviation 
of  6  percent  with  respect  to  the  modified  base  pressure  theory. 
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NOMENCLATURE 

Area  of  nozzle  throat 

Area  of  diffuser 

Area  of  nozzle  exit  (Fig.  1) 

Width  of  mixing  zone  (Fig.  8b) 

Inviscid  jet  boundary  Crocco  number,  V  1  -  Tj/Tt 

Rate  of  growth  of  mixing  zone 

Location  of  mixing  zone  relative  to  inviscid  jet 
boundary  (Fig.  12) 

Length  along  inviscid  jet  boundary  from  nozzle 
exit  to  diffuser  wall 

Mach  number 

Mach  number  of  inviscid  jet  boundary 

Static  pressure 

Base  or  cell  pressure  (Fig.  1) 

Total  pressure 

Total  pressure  downstream  of  a  normal  shock  wave 

Static  pressure  on  diffuser  wall 

Radius  of  diffuser 

Radius  of  nozzle  exit  (Fig.  1) 

Radius  of  total  pressure  probe  (Fig.  2b) 

Gas  total  temperature 

Inviscid  jet  boundary  static  temperature 

Velocity  anywhere  in  mixing  zone  (Fig.  12b) 

Velocity  at  inner  edge  of  mixing  zone  (Fig.  12b) 

Distance  from  nozzle  exit  plane  parallel  to  nozzle 
centerline  (Fig.  2) 

Distance  perpendicular  to  nozzle  centerline  (Fig.  2) 
Inviscid  jet  boundary  recompression  turning  angle 
Ratio  of  specific  heats 


ix 


A  E  DC-  T  D  R-64-3 


77  Non-dimensional  mixing  zone  coordinate  (Refs.  7, 

8,  and  9;  Fig.  12b) 

rjM  Non-dimensional  location  of  inner  edge  of  mixing 

zone  (Fig.  12b) 

7?s  Non-dimensional  location  of  stagnating  streamline 

(Fig.  12b) 

qx  Non-dimensional  location  of  inviscid  jet  boundary 

(Fig.  12b) 

0ne  Nozzle  exit  wall  angle  relative  to  nozzle  centerline 

SUBSCRIPTS 

exp  Experimental 

max  Maximum 

min  Minimum 

theo  Theoretical 
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1.0  INTRODUCTION 


The  ejector  research  program  conducted  in  the  Rocket  Test  Facility 
(RTF),  Arnold  Engineering  Development  Center  (AEDC),  Air  Force  Sys¬ 
tems  Command  (AFSC),  has  been  primarily  of  an  experimental  nature. 

The  purpose  of  the  overall  research  program  is  to  define  the  various 
factors  which  influence  the  performance  of  axisymmetric  ejector  sys¬ 
tems.  The  results  are  presented  in  Refs.  1  through  6.  Based  on  these 
experimental  results,  simple  empirical  methods  were  developed  for 
estimating  the  performance  of  axisymmetric  ejector  systems.  In  this 
report,  the  two-dimensional  base  pressure  theory  for  isoenergetic  mix¬ 
ing  and  negligible  initial  boundary  layer,  developed  by  Korst  (Ref.  7),  is 
extended  to  improve  the  applicability  of  this  theory  to  predict  the  mini¬ 
mum  cell  pressure  ratio  produced  by  axisymmetric  ejector  systems. 

The  two-dimensional  supersonic  base  pressure  theory  developed  by 
Korst  consists  of  two  independent  theories  which,  when  solved  simul¬ 
taneously,  yield  the  base  pressure.  The  criterion  for  the  simultaneous 
solution  is  that  the  total  pressure  on  the  dividing  streamline  in  the  mix¬ 
ing  zone  must  equal  the  peak  recompression  static  pressure  in  the  im¬ 
pingement  zone.  In  Ref.  7  these  theories  were  developed  for  two- 
dimensional  mixing  along  a  uniform  inviscid  jet  flow  field.  The  total 
pressure  on  the  dividing  streamline  in  the  mixing  zone  is  determined 
by  applying  overall  momentum  and  conservation  of  mass  flow  relations 
to  the  mixing  zone  using  a  velocity  distribution  derived  by  simplifying 
and  solving  the  basic  equations  for  turbulent  flow.  The  peak  recom¬ 
pression  static  pressure  is  determined  by  turning  the  inviscid  flow  field 
parallel  to  the  impingement  surface  by  means  of  a  two-dimensional  shock 
wave. 

In  Ref.  7  it  is  shown  that,  for  axisymmetric  ejector  systems,  Korst1  s 
two-dimensional  supersonic  base  pressure  theory  predicts  minimum  cell 
pressure  ratios  greater  than  experimental  values  (1.60  <  Mj  <  3.50)  with 
the  deviation  increasing  as  the  area  ratio  of  the  driving  nozzle  increases. 
The  basic  assumptions  of  the  theory  are  reviewed  (Ref.  8)  with  regard 
to  their  compatibility  with  the  conditions  that  are  produced  by  an  axisym¬ 
metric  ejector  system.  The  results  of  that  investigation  of  interest  here 
are  as  follows: 

1.  The  theoretical  total  pressure  on  the  dividing  streamline.. _ 
agrees  with  the  measured  peak  recompression  static 
pressure  if  it  is  assumed  that  the  stagnation  process 
is  isentropic  and  that  the  mixing  is  two-dimensional. 
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2.  The  theoretical  peak  recompression  static  pressure 

obtained  by  applying  two-dimensional  shock  wave  theory 
based  on  a  free-jet  Mach  number  determined  by  the  ratio 
pc/pt  is  shown  to  disagree  seriously  with  the  measured 
peak  recompression  static  pressure. 

These  conclusions  have  also  been  experimentally  verified  by  data 
obtained  at  the  AEDC.  In  addition  to  the  above-mentioned  considera¬ 
tions,  the  axisymmetric  effects  on  the  characteristics  of  a  mixing  zone 
have  been  theoretically  investigated  in  Refs.  8  and  9.  Application  of  the 
theory  presented  in  these  references  to  the  typical  ejector  systems 
studied  in  Refs.  1  through  6  shows  that  the  axisymmetric  effects  on  mix¬ 
ing  are  negligible  for  these  systems. 

The  above  discussion  leads  to  the  conclusion  that  the  major  source 
of  error  in  the  application  of  Korst's  two-dimensional  supersonic  base 
pressure  theory  to  axisymmetric  ejector  systems  is  in  the  method  used 
for  predicting  the  peak  recompression  pressure.  An  empirical  method 
for  estimating  the  peak  recompression,  based  on  relatively  limited  ex¬ 
perimental  data  (2.57  <  Mj  <  3.41),  is  presented  in  Ref.  8.  The  corrected 
base  pressure  theory  is  then  shown  to  predict  very  accurately  the  mini¬ 
mum  cell  pressure  ratio  of  a  limited  number  of  axisymmetric  ejector 
systems.  In  this  report,  a  theoretical  method  of  estimating  the  peak 
recompression  pressure  is  developed,  and  the  modified  base  pressure 
theory  is  applied  to  a  wide  range  of  axisymmetric  ejector  systems 
(1.60  <  Mj  <  7.4)  . 

The  recompression  theory  developed  in  this  report  is  applicable  to 
any  system  involving  the  impingement  of  a  free  jet  on  a  solid  surface. 


2.0  RECOMPRESSION  MECHANISM 


The  basic  premise  of  this  analysis  is  that  the  recompression  static 
pressure  distribution  is  a  function  of  the  combined  inviscid  flow  field 
and  the  mixing  zone  flow  field.  The  following  discussion  treats  the 
theoretical  inviscid  flow  field  and  the  superposition  of  a  mixing  zone. 


2.1  THEORETICAL  INVISCID  JET  FLOW  FIELD 

A  typical  axisymmetric  ejector  system  is  shown  in  Fig.  1.  The 
inviscid  flow  field  of  the  jet  emanating  from  the  nozzle  exit  can  be 
determined  by  applying  the  method  of  characteristics  solution  to  the 
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general  potential  flow  field  equations  using  the  known  nozzle  exit  flow 

conditions  and  the  cell  pressure  ratio  This  method  is  adequately- 

covered  in  Ref.  10  for  both  isentropic  (irrotational)  and  rotational  (total 
pressure  gradient)  flow.  From  the  practical  point  of  view,  the  assump¬ 
tion  of  isentropic  flow  simplifies  the  flow  field  calculation;  however,  this 
obviously  is  a  questionable  assumption  if  the  nozzle  exit  flow  is  known 
to  be  highly  rotational  because  of  the  contour  of  the  nozzle.  Even  if  the 
nozzle  flow  is  isentropic,  the  jet  flow  field  cannot  be  treated  in  general 
as  isentropic  flow  because  of  the  presence  of  a  shock  wave  referred  to 
herein  as  the  jet  boundary  shock  wave.  The  jet  boundary  shock  wave  is 
formed  by  the  coalescence  of  infinitesimal  compression  waves  which 
result  from  the  curvature  of  the  axisymmetric  constant  pressure  jet 
boundary  (Ref.  11).  The  strength  (change  in  entropy)  of  the  jet  boundary 
shock  wave  increases  with  jet  radius.  Fortunately,  for  most  ejector 
systems,  the  jet  impinges  on  the  diffuser  wall  before  the  boundary  shock 
wave  can  develop  significant  strength  and  therefore  may  be  treated  as 
isentropic.  This  is  also  indicated  by  the  fact  that  the  location  of  the 
inviscid  jet  boundary  can  be  accurately  calculated  neglecting  the  boundary 
shock  wave. 


2.2  EXPERIMENTAL  VERIFICATION  OF  THEORETICAL  INVISCID  FLOW  FIELD 

On  the  basis  of  the  previous  discussion  the  jet  flow  field  is  assumed 
to  be  isentropic  and  is  theoretically  treated  by  the  method  presented  in 
Ref.  12.  To  verify  this  theoretical  method,  lines  of  constant  Mach  num¬ 
ber  in  the  free  jets  of  two  ejector  systems  were  experimentally  deter¬ 
mined  from  total  pressure  measurements  and  the  assumption  of  isen¬ 
tropic  flow.  Thus,  the  total  pressure  anywhere  in  the  flow  field  was 
assumed  to  be  equal  to  the  nozzle  plenum  total  pressure.  Since  a  total 
pressure  probe  in  a  supersonic  stream  measures  the  total  pressure 
downstream  of  a  normal  shock  wave,  the  ratio  of  measured  total  pres¬ 
sure  to  nozzle  plenum  total  pressure  can  be  used  to  determine  the  free- 
stream  Mach  number. 

The  two  ejector  systems  tested  used  18-deg,  half-angle  conical 
driving  nozzles  having  area  ratios  of  10.  848  and  25.  000,  which  are 
described  in  Ref.  4.  The  diffuser  diameter  for  both  configurations  was 
10.  19  in.  The  total  pressures  in  the  internal  portion  of  the  jet  flow  field 
were  measured  by  a  fixed  rake  (Fig.  2a),  and  the  driving  nozzle  was 
axially  translatable  relative  to  the  fixed  rake;  thus  measurements  could 
be  obtained  at  various  axial  stations.  A  complete  description  of  the  rake 
and  movable  nozzle  system  is  presented  in  Ref.  6.  The  driving  fluid  was 
unheated  air  at  a  temperature  of  about  80 °F. 
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The  theoretical  flow  field,  based  on  the  assumption  of  isentropic 
flow,  is  shown  in  Figs.  3  and  4  to  agree  well  with  the  experimental 
results  for  the  internal  portion  of  the  free  jet.  The  theoretical  flow 
field  was  calculated  by  the  IBM  7070  computer.  The  experimental 
results  presented  in  Figs.  3  and  4  were  obtained  by  cross -plotting  the 
original  experimental  data  in  such  a  manner  as  to  obtain  lines  of  con¬ 
stant  Mach  number. 

Experimental  flow  field  measurements  in  the  vicinity  of  the  jet 
boundary  were  made  with  a  single  total  pressure  probe,  approximately 
aligned  with  the  flow  direction  as  shown  in  Fig.  2b.  The  jet  flow  field 
in  this  region  cannot  be  assumed  to  be  irrotational  (constant  total  pres¬ 
sure)  because  of  the  mixing  zone  and  the  jet  boundary  shock  wave. 
Therefore,  the  measured  total  pressure  cannot  be  used  to  determine 
Mach  number  without  making  further  assumptions  concerning  the  nature 
of  the  actual  flow  field  or  by  measuring  the  static  pressures.  Static 
pressure  measurements  in  this  portion  of  the  flow  field  are  impractic¬ 
able  by  conventional  methods  because  of  the  non-uniformity  of  the  flow 
field.  To  avoid  these  difficulties,  the  ratio  of  the  experimentally  meas¬ 
ured  total  pressure  to  the  nozzle  total  pressure  is  compared  in  Fig.  5 
with  the  theoretical  total  pressure  that  the  probe  should  sense.  This 
comparison  shows  that  the  theoretical  location  of  the  jet  boundary  shock 
wave  and  the  flow  conditions  upstream  of  the  boundary  shock  wave  agree 
well  with  the  experimental  results.  It  is  also  shown  in  Fig.  5  that  in  the 
region  between  the  jet  boundary  shock  wave  and  the  inviscid  jet  boundary, 
the  theoretical  flow  field  is  in  error  because  of  the  previously  mentioned 
presence  of  the  mixing  zone  and  the  nonisentropic  nature  of  the  jet 
boundary  shock  wave.  For  the  particular  ejector  systems  used  in  these 
experiments,  the  theoretical  maximum  total  pressure  loss  through  the 
jet  boundary  shock  wave  is  less  than  5  percent.  Therefore,  the  mixing 
zone  is  the  primary  cause  of  the  disagreement  between  the  theory  and 
experiment  in  this  region  for  these  particular  configurations. 

In  summary,  it  can  be  concluded  that  the  jet  flow  field,  out  to  the 
impingement  shock  wave,  is  accurately  predicted  by  applying  the  method 
of  characteristics  solution  to  the  general  isentropic  flow  equations  except 
in  the  region  of  the  mixing  zone. 


2.3  THEORETICAL  RECOMPRESSION  STATIC  PRESSURE  DISTRIBUTION 

The  flow  field  approaching  the  diffuser  wall  is  rotational  in  the  mix¬ 
ing  zone  portion  and  non-uniform  in  the  inviscid  portion.  The  turning 
of  such  a  flow  field  parallel  to  a  wall  can  be  treated  theoretically  (neg¬ 
lecting  boundary  layer)  to  determine  the  static  pressure  distribution 
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along  the  wall  by  considering  small  increments  in  the  approaching  flow 
to  be  uniform  at  some  average  Mach  number,  flow  direction,  and  total 
pressure  level.  The  impingement  shock  wave  can  be  treated  as  two- 
dimensional  for  most  systems  of  practical  interest  since  it  is  always  at 
a  large  radius  from  the  centerline  of  the  nozzle.  A  sketch  of  a  repre¬ 
sentative  first-order  theoretical  recompression  flow  field  is  presented 
in  Fig.  6.  This  is  referred  to  as  a  first-order  theory  since  the  weaker 
or  second-order  interaction  waves  are  neglected.  The  calculation  pro¬ 
cedure  involves  the  following: 


The  pressure  Pn  is  produced  by  the  turning  (two- 
dimensionally  of  the  average  flow  conditions  in  region 

1  through  the  angle  |^—  *  . 

The  average  flow  in  region  2  turns  through  the  angle 


(A±&)  -  s, 


in  order  to  satisfy  the  interaction 


conditions  with  the  downstream  flow  from  region  1. 

The  interaction  conditions  are:  equal  static  pressures 
(p2i  =  p22)  and  parallel  flow.  Thus,  the  angle  8,  can  be 
determined  by  a  trial  and  error  procedure.  The  magni¬ 
tude  of  the  angle  8,  is  typically  about  1  deg  or  less, 
depending  on  the  size  of  the  region.  Therefore,  the  pres¬ 
sure  P2i  is  determined  by  an  isentropic  compression  from 
Pu  through  the  angle  8,  . 


3.  The  pressure  p31  is  produced  by  an  isentropic  compression 
from  p21  through  the  angle  8,.  The  pressure  p32  is  pro¬ 
duced  by  an  isentropic  compression  from  p22  through  the 
angle  8,  .  When  second-order  compression  and  expansion 
waves  are  neglected,  p32  =  p31  • 


4.  The  pressure  p41  is  produced  by  an  isentropic  compres¬ 
sion  from  p31  through  the  angle  s2  .  The  pressure  p42  is 
produced  by  an  isentropic  compression  from  p32  through 
the  angle  S2  .  The  pressure  p43  is  produced  by  the  turning 
of  the  average  flow  in  region  3  through  the  angle  82  in  such 
a  manner  as  to  satisfy  the  previously  stated  interaction 
conditions.  When  second-order  waves  are  neglected, 

P*3  _  P<2  ~  P4I  ' 

5.  The  pressure  p  is  produced  by  an  isentropic  compression 
from  p41  through  the  angle  82  •  The  pressure  p52  is  produced 
by  an  isentropic  compression  from  p42  through  the  angle  82  . 
The  pressure  p53  is  produced  by  isentropic  compression 
from  p43  through  the  angle  82  .  When  second-order  waves 
are  neglected,  psi  =  p52  =  p53  • 
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Steps  1  through  5  treat  the  theoretical  recompression  of  the  mixing 
zone  portion  of  the  approaching  flow  field.  The  theory  shows  that  the 
static  pressure  distribution  along  the  diffuser  wall  is  of  an  increasing 
nature  from  pu  to  psl  -  produced  by  a  series  of  isentropic  compressions. 

It  is  also  shown  that  the  approaching  flow  does  not  turn  parallel  to  the 
diffuser  wall,  initially,  but  through  the  angle  (/3  -  8) .  The  flow  is  then 
turned  parallel  to  the  diffuser  wall  by  an  isentropic  compression  through 
the  angle  8  . 

The  theoretical  recompression  of  the  inviscid  portion  of  the  approach¬ 
ing  flow  (constant  total  pressure  and  increasing  Mach  number)  can  be 
treated  in  a  similar  manner. 

6.  The  pressure  p64  is  produced  by  the  turning  of  the  average 

flow  conditions  in  region  4  through  the  angle  j^3  *  —  +  S3J 

in  such  a  manner  as  to  satisfy  the  interaction  conditions 
with  the  downstream  flow  from  region  3. 

The  interaction  conditions  require  the  downstream  flow  from  region  3 
to  expand  isentropically  through  the  angle  S3  to  the  pressure  p63  because 
of  the  characteristics  of  the  two-dimensional  shock  wave  theory  shown, 
for  a  typical  case,  in  Fig.  7  (decreasing  downstream  static  pressure 
with  Mach  number).  Therefore,  the  pressure  p62  is  produced  by  an 
isentropic  expansion  from  p52  through  the  angle  <53  .  The  pressure  p61 
is  produced  by  an  isentropic  expansion  from  p51  through  the  angle  S3  . 

When  second-order  waves  are  neglected,  p64  =  p63  =  p62  =  p6i  . 

7.  The  pressure  p71  is  produced  by  an  isentropic  expansion 
from  p6i  through  the  angle  S3  .  In  a  similar  manner  the 
pressures  p72  ,  p73  ,  and  p74  are  determined.  When  second- 
order  waves  are  neglected,  p74  =  p73  =  p72  =  p71  • 

8.  The  above  procedure  is  repeated  for  regions  5,  6,  7,  etc. 

Steps  6  through  8  show  the  theoretical  diffuser  wall  static  pressure 
distribution  to  be  of  a  decreasing  nature  from  p51  .  This  analysis  also 
shows  that  the  approaching  flow  does  not  turn  parallel  to  the  diffuser 
wall,  initially,  but  through  the  angle  (yS  +  8)  .  The  flow  is  then  turned 
parallel  to  the  diffuser  wall  by  an  isentropic  expansion  through  the 
angle  8 . 

The  error  in  this  theory  caused  by  neglecting  the  second-order 
waves  can  be  estimated  in  the  following  manner.  The  flow  interaction 
requirements  which  determine  the  impingement  shock  wave  are  equal 
static  pressures  and  flow  direction.  However,  the  Mach  numbers  in 
the  two  regions  produced  by  the  interaction  will  be  different  because 
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of  the  differences  in  the  total  pressure  levels.  As  an  example,  the 
Mach  number  in  the  region  where  p2l  exists  may  be  1.  60  and  the  Mach 
number  in  the  region  where  p22  exists  may  be  2.  6,  depending  on  the 
size  of  the  flow  field  increment.  If  St  =  1.0  deg,  then  from  Fig.  8 
for  an  isentropic  compression, 


Pn 


P22 

Since  p2l  =  p22 ,  second-order  waves  are  necessary  to  satisfy  the  inter¬ 
action  conditions.  Neglecting  second-order  waves  results  in  an  error 
of  +1.  91  percent  in  p32  relative  to  p31  .  The  difference  in  Mach  number 
used  in  this  example  is  relatively  large;  a  typical  Mach  number  differ¬ 
ence  would  be  0.25,  and  the  corresponding  error  in  p32  relative  to  p31 
would  be  0.  42  percent.  The  curves  in  Fig.  8  show  the  relative  insensi¬ 
tivity  of  the  pressure  ratio  with  Mach  number  for  Mach  numbers  greater 
than  1.  1.  Thus,  the  first-order  theoretical  recompression  analysis  is 
accurate  for  at  least  that  portion  of  the  recompression  static  pressure 
distribution  from  the  stagnating  streamline  to  a  point  slightly  beyond 
the  peak  static  pressure  location. 

For  typical  flows  considered,  the  theoretical  peak  recompression 
static  pressure  is  psi  and  is  shown  in  step  5  to  be  equal  to  p53  .  The 
pressure  p53  is  determined  by  the  average  flow  conditions  in  the  vicinity 
of  the  inner  edge  of  the  mixing  zone  (region  3).  Since  the  initial  turning 
angle  of  the  mixing  zone  portion  of  the  approaching  flow  field  is  (/3  -  8) 
and  since  the  initial  turning  angle  of  the  inviscid  portion  of  the  approach¬ 
ing  flow  field  is  ,(/5  +  8)  ,  then  the  streamline  along  the  inner  edge  of  the 
mixing  zone  must  turn  through  the  angle  /3  .  Thus,  the  end  of  the  com¬ 
pression  process  and  beginning  of  the  expansion  process  is  determined 
by  the  location  of  the  inner  edge  of  the  mixing  zone.  As  a  result,  the 
peak  recompression  static  pressure  is  equal  to  the  static  pressure  down¬ 
stream  of  the  impingement  shock  wave  required  to  turn  the  flow  along 
the  inner  edge  of  the  mixing  zone  parallel  to  the  diffuser  wall. 


1.051 

1.071 


2.4  SIMPLIFIED  THEORETICAL  RECOMPRESSION  STATIC  PRESSURE  DISTRIBUTION 

The  first-order  theoretical  recompression  static  pressure  distribu¬ 
tion  involves  a  trial  and  error  procedure  to  determine  S  and  a  scale 
drawing  of  the  flow  field  to  determine  the  shock  wave  arrangement.  The 
recompression  static  pressure  distribution  can  be  approximated  by  ignor¬ 
ing  the  actual  impingement  shock  wave  configuration  and  the  flow  field 
interactions.  The  procedure  is  to  determine  the  flow  conditions -approach¬ 
ing  the  diffuser  wall  at  each  point  and  to  turn  each  streamline  parallel  to 
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the  wall  through  an  oblique  shock.  Two  types  of  error  exist  in  this 
simplified  analysis.  One  error  is  associated  with  the  orientation 
of  the  flow  and  the  other  is  associated  with  the  magnitude  of  the 
static  pressure.  The  orientation  error  is  not  significant  if  the  angle, 
relative  to  the  diffuser  wall,  of  the  approaching  flow  field  is  approxi¬ 
mately  equal  to  the  Mach  angle  of  the  flow  downstream  of  the  impinge¬ 
ment  shock  wave.  This  is  usually  the  case  for  many  ejector  systems. 
The  error  in  the  magnitude  of  the  static  pressure  can  be  estimated  in 
the  following  manner.  Consider  the  streamline  approaching  the  diffuser 
wall  that  has  the  same  flow  conditions  as  the  average  flow  in  region  2 
(Fig.  6).  Let  p3  be  the  static  pressure  downstream  of  the  two- 
dimensional  shock  wave  required  to  bend  this  streamline  parallel  to  the 

wall.  The  turning  angle  would  be  (  —  *  .  The  static  pressure  p3 

can  be  approximated  by  considering  that  this  streamline  turns  by  means 

of  a  two-dimensional  shock  wave  through  the  angle  —  S3J  and 

then  isentropically  compresses  through  the  angle  to  the  static  pres¬ 
sure  p3 .  The  error  in  p3  relative  to  p3  for  =  1  deg  is  shown  to  be 
very  small  in  Fig.  9  for  typical  flow  conditions.  In  steps  2  and  3 
(section  2.  3)  the  static  pressure  p32  is  shown  to  be  determined  by  the 
same  process  as  that  used  to  determine  p3  .  Therefore, 

P3  =  Pj  =  p32  =  p31  •  Thus,  the  error  in  the  magnitude  of  the  static  pres¬ 
sure  obtained  by  this  simplified  analysis  is  on  the  order  of  that  produced 
by  neglecting  second-order  waves  (section  2.  3). 


2.5  SUPERPOSITION  OF  MIXING  ZONE  ON  INVISCID  FLOW  FIELD 

The  theoretical  recompression  of  a  known  approaching  flow  field 
is  presented  in  sections  2.  3  and  2.4.  However,  for  axisymmetric  jets, 
the  mixing  takes  place  in  a  region  of  non-uniform  flow  both  alone  and 
across  the  mixing  zone.  The  greatest  flow  field  variation  is  in  the 
direction  across  the  mixing  zone,  as  shown  in  Figs.  10a  and  11a. 

These  flow  conditions  violate  the  basic  assumptions  of  uniform  flow  and 
constant  static  pressure  made  in  Ref.  7  to  theoretically  determine  the 
velocity  distribution  in  the  mixing  zone.  Therefore,  a  direct  super¬ 
position  of  this  theoretical  mixing  zone  velocity  profile  on  the  theoretical 
inviscid  flow  field  of  an  axisymmetric  jet  is  unjustifiable.  However,  as 
previously  stated  in  this  report,  only  the  peak  recompression  static  pres 
sure  is  of  importance  in  Korst's  supersonic  base  pressure  theory.  In 
section  2.  3  it  is  theoretically  shown  that  the  peak  recompression  static 
pressure  is  produced  by  the  turning,  parallel  to  the  diffuser  wall,  of  the 
flow  field  along  the  inner  edge  of  the  mixing  zone. 

A  comparison  of  experimental  recompression  static  pressure 
distribution  with  the  simplified  theoretical  distribution  based  on  the 
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isentropic  inviscid  flow  field  is  presented  in  Figs.  10b  and  lib.  The 
theoretical  distribution  is  shown  to  be  in  error  due  to  neglecting  the 
mixing  process.  However,  the  theoretical  distribution  is  shown  to 
agree  reasonably  well  with  the  experimental  data  in  the  region  down¬ 
stream  of  the  apparent  mixing  zone,  thus  verifying  the  simplified 
theoretical  analysis. 


The  region  of  maximum  influence  of  the  mixing  process  relative 
to  the  inviscid  flow  field  can  be  estimated  from  an  idealized  concep¬ 
tion  of  the  recompression  flow  phenomena.  Based  on  the  results  pre¬ 
sented  in  Figs.  10b  and  lib,  an  idealized  recompression  flow  model 
was  developed  as  shown  in  Figs.  12a  and  b.  The  region  of  influence 
of  the  mixing  zone  on  the  recompression  static  pressure  distribution 
can  be  determined  in  the  following  manner: 

a.  Impingement  location  of  the  inner  edge  of  the  mixing  zone 


For  turbulent  mixing 
b 

-  =  C 

rn  e 

Substituting  Eq.  (2)  into  Eq.  (1)  yields 
=  kc  /  \ 

rne  sin  fi  V  rne  / 


(1) 

(2) 

(3) 


Values  for  k  and  c  (from  Abramovich,  Ref.  13)  are  derived  in 

£  * 

Ref.  9  and  presented  in  Fig.  13.  Values  for  y-jy  and  /3  are 
obtained  from  the  theoretical  inviscid  jet  boundary. 


b.  Impingement  location  of  stagnating  streamline 


The  thickness  of  the  mixing  zone  at  the  impingement  point 
of  the  stagnating  streamline  is 


(4) 


by  geometry 


Ax2 


r 


n  e 


by  definition  (Ref.  9) 


27?m 


(5) 

(6) 


9 


A  E  DC-T  D  R-64-3 


Substituting  Eqs.  (4)  and  (6)  into  Eq.  (5)  and  solving 
for  f*2-  yields 


k  c 


/  ij  \(ru  -v,  \ 

{  r„  eJ\VM~Voa) 


1  +  ^  (^^~) 

tzn/3  \  T)m 


Location  of  outer  edge  of  mixing  zone 

The  equation  for  the  thickness  of  the  mixing  zone  is 


(&) 


2  '  cos  /3  - 


Axx 


Ax 

By  geometry,  —  1  can  be  approximated  by  the  following 
equation 


(7) 


(8) 


Axj 


(1  -  k) 

(— )  - 

(  AX2)  sin/sl 

_ 

V  rne  / 

\  rne  / 

tan 


P 


O) 


Substituting  Eq.  (8)  into  Eq.  (9)  and  solving  for  ~~l  ■  yields 

Ax, 


(  i  -  k) 

C  | 

(*  i  ) 

l  rn  e  j 

1  —  [  (  1  —  k  )  c 

+  tan/3]  (4~)  cos/3 

(l 

-k) 

p  ~ 

c  +  tan  —  J 

(10) 


The  location  tjs  of  the  stagnating  streamline  is  determined  by 
the  condition  that  the  total  pressure  of  the  stagnating  stream¬ 
line  equals  the  peak  recompression  static  pressure. 

The  theoretical  mixing  zone  width  determined  by  Eqs.  (3),  (7),  and 
(10)  is  shown  in  Figs.  10b  and  lib  to  agree  very  well  with  the  apparent 
experimental  width  of  the  mixing  zone.  It  is  also  shown  in  Figs.  10b 
and  lib  that  the  theoretical  recompression  static  pressure  determined 
by  the  inviscid  flow  conditions  along  the  inner  edge  of  the  mixing  zone 
agrees  in  magnitude  with  the  experimental  peak  recompression  static 
pressure. 


2.6  EXPERIMENTAL  VERIFICATION  OF  THE  THEORETICAL  PEAK 
RECOMPRESSION  STATIC  PRESSURE 

In  section  2.  3,  it  was  theoretically  shown  that  the  peak  recompres¬ 
sion  static  pressure  is  determined  by  turning,  parallel  to  the  diffuser 
wall,  the  inviscid  flow  along  the  inner  edge  of  the  mixing  zone.  Theo¬ 
retical  peak  recompression  static  pressures  are  compared  with  experi¬ 
mental  values  in  Figs.  14  and  15.  The  theoretical  values  were  obtained 
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by  using  the  method  of  characteristics  solution  to  the  general  potential 
isentropic  flow  equations  to  determine  the  inviscid  flow  field,  based  on 

the  experimental  (-jp-)  and  Eq.  (3)  to  determine  the  location  of  the  inner 


edge  of  the  mixing  zone.  The  experimental  data  were  obtained  from 
three  ejector  systems  geometrically  described  in  each  figure.  The 


minimum  cell  pressure  ratio 


for  each  configuration  was  varied 


by  varying  the  total  pressure  (pt)  level  of  the  driving  nozzle.  The 
resulting  variation  of  (pc/pt)  with  pt  is  shown  in  Fig.  16  and  is  due  to  a 
Reynolds  number  effect  described  in  Ref.  4.  The  Reynolds  number 
determines  the  character  of  the  nozzle  boundary  layer  which  then  in¬ 
fluences  the  free-jet  mixing  process,  thus  producing  a  new  equilibrium 
cell  pressure.  The  nozzle  boundary  layer  influences  the  jet  mixing 
process  in  two  ways. 


1.  The  velocity  profile  in  the  mixing  zone  varies  in  shape 
with  distance  from  the  nozzle  exit.  However,  the 
asymptotic  profile  is  the  same  as  the  undisturbed  mix¬ 
ing  profile. 

2.  The  mixing  zone  has  a  finite  thickness  at  the  nozzle 
exit  and,  consequently,  is  thicker,  at  any  point,  than 
the  undisturbed  mixing  zone. 


Theoretically,  the  two  characteristics  of  the  mixing  zone,  which  are 
involved  with  the  peak  recompression  static  pressure,  are  k  and  b.  The 
relative  location  (k)  is  a  function  of  the  velocity  profile  in  the  mixing  zone. 
The  width  (b)  of  the  mixing  zone  at  the  point  of  impingement  is  a  func¬ 
tion  of  the  thickness  of  the  mixing  zone  at  the  nozzle  exit.  For  most 
ejector  systems  of  practical  interest,  the  length  of  the  mixing  zone  is 
sufficient  to  allow  the  velocity  profile  to  transform  into  the  undisturbed 
profile.  Thus,  the  relative  location  (k)  is  the  same  as  for  the  undis¬ 
turbed  case.  However,  the  width  of  the  mixing  zone  (b)  for  most 
ejector  systems  of  practical  interest  will  be  thicker  than  the  undisturbed 
case  and  will  produce  a  lower  peak  recompression  static  pressure  than 
the  undisturbed  case.  This  is  experimentally  verified  in  Figs.  15b  and 
c  by  comparing  the  theoretical  peak  recompression  static  pressure  for 
an  undisturbed  mixing  zone  with  the  experimental  values  obtained  at  low 
total  pressure  levels  (low  Reynolds  number).  The  comparison  in  Fig.  15 
shows  that  the  theoretical  peak  recompression  static  pressure  deter¬ 
mined  by  the  flow  conditions  along  the  inner  edge  of  the  mixing  zone 
agrees  to  within  10  percent  of  the  experimental  values  of  the  maximum 
total  pressure  level.  Included  in  Figs.  14  and  15  is  the  peak  recompres¬ 
sion  static  pressure  that  would  exist  if  there  were  no  mixing.  This  is  the 
recompression  static  pressure  produced  by  the  turning  of  the  inviscid  jet 
boundary  streamline  parallel  to  the  diffuser  wall. 
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Since  the  major  factor  which  determines  the  peak  recompression 
static  pressure  is(Pc/Pt),  theoretical  peak  recompression  static  pres¬ 
sures  are  compared  with  experimental  values  in  Fig.  14  as  a  function 
of  (Pc/Pt).  These  are  the  same  experimental  values  as  presented  in 
Fig.  15.  Included  in  Fig.  14c  are  experimental  peak  recompression 
static  pressures  obtained  with  a  small  amount  of  bleed  flow  into  the 
cell  region.  The  bleed  flow  caused  the  cell  pressure  ratio  to  increase 
to  a  new  equilibrium  value  and  a  new  peak  recompression  static  pres¬ 
sure.  The  theoretical  peak  recompression  static  pressures  determined 
by  the  flow  conditions  along  the  inner  edge  of  the  mixing  zone  are  shown 
in  Fig.  14c  to  agree  very  well  with  the  new  experimental  values.  Fur¬ 
ther  experimental  verification  of  the  recompression  theory  is  presented 
in  Table  1,  in  which  the  theory  is  compared  with  data  obtained  from 
Ref.  8. 

From  the  previous  discussion  it  can  be  concluded  that  the  peak  re¬ 
compression  static  pressure  determined  by  the  flow  conditions  along  the 
inner  edge  of  the  mixing  zone  agrees  with  the  experimental  results  to 
within  ±10  percent  for  those  cases  in  which  the  nozzle  exit  boundary  layer 
can  be  neglected. 


3.0  THEORETICAL  MINIMUM  CELL  PRESSURE  RATIO 
OF  AX  SYMMETRIC  EJECTOR  SYSTEMS 


Korst  (Ref.  7)  has  developed  the  two-dimensional,  supersonic  base 
pressure  theory  for  negligible  initial  boundary  layer.  When  trying  to 
apply  this  theory  to  an  axisymmetric  ejector  system,  most  of  the  basic 
assumptions  (discussed  in  Ref.  8)  are  clearly  violated.  Also  in  Refs.  8 
and  9  it  is  shown  that  the  total  pressure  on  the  dividing  streamline  is, 
for  most  practical  ejector  systems,  independent  of  axisymmetric  effects. 
The  non -uniformity  of  the  inviscid  flow  would  be  expected  to  produce  a 
significant  deviation  from  the  two-dimensional  value  of  the  total  pressure 
on  the  dividing  streamline;  however,  in  Ref.  8  and  indirectly  by  the  data 
in  this  report  no  appreciable  deviation  was  noted.  A  possible  explana¬ 
tion  of  this  is  the  following.  The  non-uniformity  of  the  inviscid  flow 
field  distorts  the  mixing  zone  in  such  a  manner  as  to  cause  the  total 
pressure  on  the  dividing  streamline  to  be  greater  than  the  two- 
dimensional  value  (Ref.  8).  The  dividing  streamline  may  then  experi¬ 
ence  a  normal  shock  before  stagnating  on  the  diffuser  wall  (Ref.  17), 
thus  reducing  the  total  pressure.  These  two  processes  tend  to  cancel 
each  other,  resulting  in  a  total  pressure  level  close  to  the  two-dimensional 
isentropic  recompression  pressure. 
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As  previously  stated,  the  major  source  of  error  associated  with 
the  application  of  the  two-dimensional  base  pressure  theory  to  an 
axisymmetric  ejector  system  is  the  assumption  of  uniform  inviscid 
flow  in  the  determination  of  the  peak  static  pressure  in  the  impinge¬ 
ment  zone.  This  error  inKorst's  base  pressure  theory  has  been 
eliminated  by  using  the  'recompression  theory  developed  herein.  The 
resulting  base  pressure  theory  is  referred  to  as  the  "modified"  base 
pressure  theory. 

A  comparison  of  theoretical  minimum  cell  pressure  ratios  with 
experimental  data  is  presented  in  Fig.  16  for  various  axisymmetric 
ejector  systems.  The  theoretical  values  were  obtained  by  using  the 
base  pressure  theory  (isoenergetic,  neglecting  initial  boundary  layer) 
presented  in  Ref.  7  modified  by  the  theory  developed  in  this  report  for 
estimating  the  peak  reconjpression  static  pressure.  Included  in 
Figs.  16a  and  b  are  theoretical  minimum  cell  pressure  ratios  using 
the  unmodified  two-dimensional  theory.  As  shown  in  Figs.  16a  through 
e,  the  modified  theory  agrees  very  well  with  experimental  results, 
except  for  the  two  ejector  configurations  in  Figs.  16d  and  e  which  were 
using  "bell"  type  nozzles.  Theoretically,  these  two  ejector  configura¬ 
tions  were  treated  as  having  isentropic  nozzles  with  uniform  exit  flow 
conditions;  this,  of  course,  is  not  correct  since  the  total  pressure  and 
velocity  vary  non-uniformly  across  the  exit.  This  is  believed  to  be  the 
reason  for  the  major  portion  of  the  difference  between  the  theoretical 
and  experimental  values  obtained  for  these  two  configurations. 

All  of  the  data  presented  in  Fig.  16  were  obtained  from  ejector  sys¬ 
tems  using  air  (y  =  1.4)  as  the  driving  fluid  and  having  nozzle  exit  flow 
angles  of  0,  15,  and  18  deg.  In  Table  2  theoretical  minimum  cell  pres¬ 
sure  ratios  are  compared  with  experimental  values  (Ref.  14)  for  an 
ejector  system  using  a  conical  nozzle  having  a  half  angle  of  7.  58  deg. 
This  system  was  operated  using  both  air  (y  =  1.40)  and  helium  (y  =  1.66) 
as  the  driving  fluids;  the  theoretical  minimum  cell  pressure  ratios  are 
shown  to  agree  very  well  with  the  experimental  values. 

In  summary,  the  theoretical  results  were  compared  with  experi¬ 
mental  values  obtained  from  34  axisymmetric  ejector  systems,  neg¬ 
lecting  "bell"  ejector  systems,  having  the  following  characteristics: 

y  =  1.4  and  1.66 

1.00  <  -Af.  <  25.00 
A* 

1.46  <  Alf-  <  150.36 
A* 

0  deg  <  dae  <  18  deg 
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The  use  of  the  modified  theory  resulted  in  minimum  cell  pressure  ratio 
values  which  deviated  from  the  experimental  values  by  a  maximum  of 
-11.8  to  +25.  2  percent,  with  a  standard  deviation  of  5,  5  percent 
(Ref.  15). 

The  modified  base  pressure  theory  also  applies  directly  to  the 
two-dimensional  case  with  a  curved  jet  boundary  as  compared  with  the 
simpler  theory  derived  by  Korst  for  straight  jet  boundaries.  The  modi¬ 
fied  base  pressure  theory  shows,  in  general,  that  the  minimum  cell 
pressure  ratio  produced  by  an  axisymmetric  ejector  system  is  a  func¬ 
tion  of  the  thickness  of  the  mixing  zone  at  the  free-jet  boundary  impinge¬ 
ment  point,  unlike  the  two-dimensional  back  step  case  which  has  a 
uniform  inviscid  flow  field. 


4,0  CONCLUSIONS 


Based  on  the  results  of  this  study,  the  following  conclusions  are 
reached: 

1.  The  peak  recompression  static  pressure  of  an  axisym¬ 
metric  jet  is  approximately  determined  by  the  inviscid 
flow  conditions  along  the  inner  edge  of  the  mixing  zone. 

2.  The  peak  recompression  static  pressure  and  therefore 
the  minimum  cell  pressure  ratio  of  an  axisymmetric 
ejector  system  is,  in  general,  a  function  of  the  thick¬ 
ness  of  the  mixing  zone  at  the  free-jet  boundary  impinge¬ 
ment  point. 

3.  The  minimum  cell  pressure  ratio  of  the  axisymmetric 
ejector  systems  investigated  can  be  estimated  to  within 
6  percent  (standard  deviation)  through  an  application 

of  Korst' s  isoenergetic,  supersonic  base  pressure  theory, 
modified  by  the  recompression  theory  presented  in  this 
report.  These  results  were  verified  for  jet  boundary 
Mach  number’s  from  1.6  to  7.  4. 
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TABLE  1 

COMPARISON  OF  THEORETICAL  WITH  EXPERIMENTAL  PEAK  RECOMPRESSION  STATIC  PRESSURE 

y  =  1.4 


An  e 

r'D 

^ne 

deg 

Nozzle 

Type 

Pc 

/  ^wmax| 

Error, 

percent 

Remarks 

A* 

rn  e 

Pt 

\  Pc  ) 

\  P°  /theo 

1.  0 

2.  00 

0 

Isentropic 

0. 05283 

3.  06 

2.  77 

-9.  48 

Data  from  Ref.  8 

1. 1065 

1.  843 

0 

M 

0. 04149 

3.  18 

3.  56 

+  10.  67 

1  T 

1.  568 

1.550 

0 

I ! 

0. 02274 

3.  68 

4.  05 

+  10.  05 

l  ! 

1.  8562 

1.440 

0 

_ i 

!  I 

0. 01501 

4.  45 

4.  86 

+9.  21 

1! 

TABLE  2 

COMPARISON  OF  THEORETICAL  WITH  EXPERIMENTAL 
MINIMUM  CELL  PRESSURE  RATIOS  OF  AXISYMMETRIC  EJECTORS 


Ad 

An  e 

$ne, 

y 

Error, 

Remarks 

A* 

A* 

deg 

V  Pt,Axp 

VPt/theo 

percent 

22.  25 

10.  71 

7.  58 

1.4 

0. 0011 

0.  000985 

-10.46 

Data 

22.  25 

10.71 

7.  58 

66(He) 

0. 00068 

0.  000712 

+4.  71 

from 
Ref.  14 

73.  06 

17.  51 

18 

1.40 

0. 00045 

0.  00045026 

+0. 0578 
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Fig.  2  Sketch  of  Experimental  Configurations 


20 


iie  Pc  *~3 

- 5  =  10.848  -=  =  65.25  —  =  0.5331  x  10  V  =  1.40 

A*  A*  Pt 


0  1.0  2.0  3.0  4.0 

(X/r  ) 
ne 

Fig.  3  Comparison  of  Theoretica I  with  Experimental  Lines  of  Constant  Mach  Number  in  a  Free  Jet,  Ane/A*  =  10.848 


A  E  DC-T  DR-64-3 


Fig.  4  Comparison  of  Theoretical  with  Experimental  Lines  of  Constant  Mach  Number  in  a  Free  Jet,  Ane/A 
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Fig.  6  Theoretical  Recompression  Flow  Field  (First  Order) 
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b.  Data  from  Ref,  16 
Fig.  16  Continued 
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